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We show that giant quasi-bound diatomic complexes, 
whose size is typically hundreds of nm, can be formed by 
intra-cavity cold diatom photoassociation or photodissocia- 
tion in the strong atom-cavity coupling regime. 
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Cold atoms exchanging single photons with cavity 
fields have been predicted to give rise to a variety of 
fascinating motional effects, both in the "good-cavity" 
(strong-coupling) regime of Rabi oscillations Q] and in 
the "bad-cavity" regime of nearly exponential decay J2|. 
Here we pose a question that has not been considered 
thus far: what happens when two identical cold atoms ex- 
change a photon in the strong-coupling ("good-cavity") 
regime during a collision or diatomic dissociation? We 
show here that this regime gives rise to a novel, hitherto 
unexplored, interplay of molecular dynamics and cavity 
QED effects. The resulting two-atom dynamics is dras- 
tically modified at interatomic separations typically ex- 
ceeding by 2 orders of magnitude those of currently in- 
vestigated cold-atom collisional resonances in magneto- 
optical traps p). The predicted modification is due to 
the possibility of diatomic quasibinding (scattering reso- 
nance) in a potential well formed by the competing ef- 
fects of the resonant dipole-dipole (RDDI) interaction 
and strong atom-cavity coupling. These effects have pre- 
viously been shown by our group to suppress energy ex- 
change between static atoms in a high-Q cavity B . The 
regime considered here drastically differs from that of 
atoms coupled by the RDDI in a " bad cavity" M or from 
that of coherently driven atoms in a cavity H . The qua- 
sibinding scattering potential formed by the RDDI and 
strong atom-cavity coupling is shown to support many 
vibrational diatom states. These states can be revealed 
in high-Q cavities by unusual spectral structures in reso- 
nance fluorescence and Raman scattering from colliding 
or dissociating pairs, as well as by sharp variation of the 
scattering cross-section as a function of their center-of- 
mass kinetic energy. 

The model which we consider is that of a pair of cold 
two-level atoms interacting with a single-mode cavity 
field. Cavity losses (damping) will be accounted for later 
on. In the center-of-mass (COM) frame, the radial part 
of the interaction Hamiltonian H can be split into two 
parts. The kinetic part is 
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is the separation between atoms A and B and I rep- 
resents the angular momentum quantum number. The 
other part is the adiabatic Hamiltonian J|] 
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Here uj c is the frequency of the cavity-mode, loa — ojb are 
the atomic transition frequencies, and ka,b are the cou- 
pling parameters of the atom A(B) to the cavity mode; 
C3 is the dipole-dipole coupling coefficient, a± and <j z 
are the atomic pseudo-spin operators and a(a') is anni- 
hilation(creation) operator of the field mode. 

The adiabatic Hamiltonian in Eq.(2) is exactly solv- 
able Q in the basis of three coupled atom-field states 
I e^,ffs,0 >, I e B ,g A ,0 > and | g A ,gB, 1 > where 
e j(dj) represents the excited (ground) state of the jth 
atom (j=A,B) while and 1 are the number of pho- 
tons in the cavity mode. We can write H a( iXi = Tu*>iXi 
(i = 1,2,3), with wi < uj2 < 0J3. As discussed in Ref. 
j§, for R -> 0, oj 3 -> w s = U A + C3/R 3 , oj 2 -» u) c , 
and w\ — ► uj a = u>a — C3/R 3 . This means that the near- 
resonant cavity mode is then decoupled from the symmet- 
ric and antisymmetric states | s, a) = 2~ 1 / 2 (| eA,gB,0 > 
± I es,gA,0 >). The symmetric-state and cavity- mode 
contributions become increasingly hybridized as R in- 
creases H: Wl ~ UJ a , LO2 — |(Ws + U! c — O), OJ3 ~ 

■n{u) s + w c + ft); with SI = y/2(nA + hb) 2 + (w» — w c ) 2 . 
The potentials oji and 0J2 exhibit a pseudo-crossing at a 
relatively large interatomic separation (~ 1000a , where 
ao is the Bohr radius) depending on the coupling strength 
of the atoms with the cavity field. The potential 0J2 has 
a minimum at the pseudo-crossing point (R c ), and so it 
can behave as a potential well supporting quasi-bound 
states. We illustrate these long-range adiabatic poten- 
tials in Fig.l for a pair of Cs atoms with levels 65*1/2 
and 6P3/2 (uja = ojb = 3.5172 x 10 14 5ec _1 ), with atom- 
cavity parameters as in the experiment of Ref. J/J and 
realistic dipole-dipole jq] coupling parameter. The larger 
the atom-field coupling strength, the deeper is the poten- 
tial well (C1J2) with the minimum point R c shifted towards 
shorter separations (Inset (a) to Fig.l). For the cavity pa- 
rameters of Ref. pj, the depth of the binding potential 
should lie in the radio-frequency regime. An increase of 
the atom-cavity coupling strength by 2 orders of mag- 
nitude will enable the binding potential u>2 to support 



bound states separated by microwave frequencies. In- 
set(b) to Fig.l shows the effect of cavity-atom detuning 
on the value of the potential u>2 at R c . For both red 
and blue detunings 102 acts as a binding potential, but 
increasing red detuning leads to more attractive u>2 and 
reduced R Cl while the opposite is true for blue detuning. 

The potentials drawn in Fig.l must be correlated to 
the short-range potentials |8|] in which the interatomic 
repulsion due to the overlap of electron charge distribu- 
tions is much larger than the atom-cavity coupling. For 
the Cs(6P 3 /2)-Cs(6Pi/ 2 ) system, the short-range limit 
(R < 28a ) of potential w\ should be correlated to the 
A 1 Y,+ potential, which at long ranges corresponds to 
the S — P3/2 asymptote ||. We interpolate the curve 
w\ so that it can merge with the model A £+ potential 
at separations comparable to the equlibrium position at 
R ~ lOao- Since the field-atom coupling strength «4/m 
varies sinusoidally along the standing-wave cavity mode, 
a pair of ultracold atoms is more likely to become quasi- 
bound by the binding potential L02 in the vicinity of an 
antinode than elsewhere. As we show in Fig.l, the in- 
ternuclear separation R for this quasi-bound state and a 
cavity mode of wavelength A c typically satisfies the crite- 
rion R << X c (the pseudo crossing separation R c ~ 10% 
of A c for Cs-Cs). Thus a high-Q cavity loaded with ul- 
tracold atoms may give rise to periodic distribution of 
giant quasibound diatoms near the cavity mode antin- 
odes. Within the range of separations R for which the 
adiabatic potential is significant, the atom-field coupling 
parameters ka and kb remain nearly constant. 

In order to elucidate the formation of quasibound 
states of two atoms by cold atomic photoassociation || 
or diatomic photodissociation inside a high-Q cavity, we 
adopt multi-channel scattering theory |l(| . Let the eigen- 
function of the Hamiltonian H = Hk + H a d with angular 
momentum I be given by ipi = ^\ = \(4>iXi)i- From the 
Schroedinger equation Htpi = Eipi, we obtain, after some 
algebra, 
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where $ = (fa fa fa) T is the solution vector and U = 
ft, diag.(wi, W2,W3) is the diagonal eigenvalue matrix, 
the nonadiabatic-coupling matrix r l > has the elements 
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higher order terms r, 

dropped, consistently with the low-energy regime. For 
the sake of convenience, Eq.(3) can be transformed into 
the following form p0[ 
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where <& 



^[TUT- 1 - ft, diag.{wi(i? 



T$, and W - . „ 
oo),u)2(R — > oo),u>3(R — > 00)}]; T is a transformation 
matrix T = exp(— J R T^'dR), with the lower limit of 
the integration at an arbitrary large separation denoted 



by Rao where all the potentials are flat; P is a 3 x 3 
diagonal matrix whose diagonal elements are the Pi = 
TrVM-E ~ Tuji(R -► 00)}, i = 1, 2, 3. 

Equation(4) represents a compact form of the scat- 
tering equations of three coupled channels, which are 
to be solved under the following boundary conditions: 
&R - 0) ~ and $a(R -v 00) ~ p- l Sin(PiR - 
l f)d mi + P l P m 1 Cos(P m R - l f)K l mi , where the reaction 
matrix K is related to the 5-matrix (see below) by 
K l = i(l - S^)(l + 5W)- 1 . Here i labels the incident 
channel, along with the angular momentum I. From now 
on we focus on s-wave (I = 0) scattering, which is ap- 
preciable for ultra-cold atoms |0, and omit the label /. 
The regular wavefunction which is a solution of the mul- 
tichannel scattering equations can then be expressed in 
matrix form flioi as 



$(i?) = J 



dR'G(R,R')W(R')<i>(R') 
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where $ is a 3 x 3 matrix with its zth column being the 
regular part of the solution vector $ for the ith incident 
channel, J and G are 3x3 diagonal matrices with the 
nth diagonal element being J n = Sin(P n R)/P n and 



G n = P T - 1 Sin(P n R)Cos(P n R') 
- p- 1 Cos(P n R)Sin(P n R') 
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These solutions must be supplemented by the analysis of 
the poles of the S-matrix 



S = piF(-P)F-\P)p-i 
which is expressed in terms of the Jost function p 
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In the absence of cavity dissipation, the bound states 
are given by the real poles of the S matrix whose energy 
E is below the lowest threshold, i.e., below huj\(R — > 00), 
so that all the corresponding channel momenta hPi(i = 
1,2,3) have positive imaginary values. The scattering 
resonances are given by the complex poles of the S ma- 
trix, with energies E above the lowest threshold, the 
imaginary part of the complex energy being negative. 
Resonances should be observed as a signature of the 
cavity-confined quasi-bound states. In Fig. 2 we show 
the variation of scattering cross-section (an) against en- 
ergy with only one (lowest) open channel, the remain- 
ing channels being closed. This variation reveals the 
cavity-induced resonances. They have an effective width 
T e ff = Tr + r c , where Tr is the resonance width in an 
ideal (lossless) cavity and T c is the rate of cavity dissipa- 
tion (cavity linewidth). For the case of Fig.l our estimate 
yields T R - 2 MHz. If we take @ k a ~ 120 MHz and 
r c ~ 40 MHz, the resonances associated with individual 
vibrational states of w 2 cannot be resolved. However, for 



cold Rydberg atoms p2| in a high-Q microwave cavity 
||13| . these resonances are shown to be resolvable due to 
the much larger ratio of ka/F c . 

One may observe the predicted resonances, e.g., by 
populating loosely bound states of the Cs2 potential A 
X E+ H that are below the dissociation threshold by an 
amount ranging from to 100 MHz, via photoexcitation 
of the ground molecular potential X 1 S+ Nonadiabatic 
effects at separations of a few nm will then produce a 
mixure of : E+ and 3 II tl states fl4|| . The alternative is 
photoassociation 0: As the two ground-state (S-state) 
cold atoms approach each other inside the cavity, they 
should be excited by a single photon. By tuning the 
frequency of this photon, the two atoms can be raised 
to the desired adiabatic potential io\ or u>2, their relative 
velocity can be controlled and determine the scattering- 
channel inputs in Eq.(5). 

In order to evaluate the characteristics of resonance flu- 
orescence associated with the discussed effect, we have 
approximated the potential ioi given in Fig.l by the 
Morse model, for which we can determine the vibrational 
states supported by that potential. The vibrational fre- 
quency u>2 V in the quasibinding potential L02 depends on 
the orientation of the molecular axis relative to the cav- 
ity axis, since the adiabatic Hamiltonian given in Eq.(2) 
is a function of the angle between the atomic dipole and 
the cavity field (via ka,b)- By applying the Landau- 
Zenner-Stueckelberg (LZS) [|10| semiclassical treatment 
of nonadiabatic coupling near the pseudocrossing of wi 
and u! 2 , we have calculated the energy-shifts of each of 
those vibrational states. The Franck-Condon spectrum 
of fluorescence for the transition from vth level of LO2 to 
the ground state (gg) of the diatom has been calculated 
for an exciting photon polarized along the X-axis, and an 
emitted photon along the Z axis, which is the direction 
of the cavity field. It can then be expressed as 



recoil energy for the transition P 3 / 2 — > Si/2- We display 
the spectrum [Eq.(9)] in Fig. 3. 

Spectral peaks in the range of the potential depth 
u>2 are indications of the cavity-confined quasi-bound di- 
atomic states. Our estimates show that for the parame- 
ters of Fig. 1 and Rcf . , the vibrational-line spacing is 
well below the cavity linewidth (r c ~ AQMHz). In order 
to resolve these lines properly, the ka/^c ratio should be 
further increas ed, e.g ., by using cold Rydberg atoms in a 
high-Q cavity yj]l3|] (see Fig. 2b). Yet, even in the ab- 
sence of optimal resolution, the spectral signature of the 
quasibinding effect in Fig. 3 is unambiguous, especially 
for the II symmetry. 

To conclude, we have demonstrated several unusual dy- 
namical features obtainable in the strong-coupling cavity 
regime of diatomic collisions or dissociation. These fea- 
tures are associated with vibrational quasibound states, 
which can be excited by intracavity photoassociation of 
colliding cold atoms via single photon absorption at a 
frequency below the normal dissociation threshold. An 
alternative is intracavity excitation of a cold dimer by 
single-photon absorption at a frequency just below the 
dissociation threshold. The prediction of cavity-induced 
giant quasi-bound cold diatomic complexes is the first 
indication that modifications of molecular collision dy- 
namics by cavity QED may be used to 'engineer' novel 
metastable states of matter, which will live as long as the 
cavity holds the photon inside. 
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with numerical computations. 
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Here C is a constant, ^ e ff is the effective linewidth as 
before, 0J2v(9) is corrected for LZS pseudocrossing effects 
and We is a weight factor, both depending on the ran- 
dom angle (9) between the Z-axis and the molecular axis: 
W e = Sin 2 (9) and We = Cos 2 (9) (IJ] for the S and 
n symmetry, respectively. The wavefunctions <f>2 V and 
4>gg represent the trth vibrational state in the uj 2 po- 
tential and the radial ground state, respectively. For 
a 4>2v state with E symmetry, the dipole selection rule 
forces 4> gg to have the angular momentum I = 1 in the 
S-S asymptotic ground state potential (-1/R 6 ). We have 
tfigg ~ Rji(y / 2[iE gg R/h) [|ll] , where j\ is the I — 1 spher- 
ical Bessel function and E gg is the ground-state (S-S) en- 
ergy, which is assumed to be of the order of one-photon 
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FIG.l Adiabatic potentials (in MHz) as a function of 
interatomic separation (in Bohr radii a^) for Cs atoms, 
u c — loa = 1.0MHz, ka — 120MHz, at short separations 
and near the pscudocrossing (equilibrium) position R c ~ 
2000ao, kb = 0.8ka- Inset(a): Depth of u>2 (in MHz) and 
position R c (in Bohr radius) as a function of atom-field 
coupling strength ka (in MHz). Inset (b): Minimum of 
<x>2 (in MHz) at R c as a function of cavity-atom detuning 
(in MHz) for red (u> c < lja) and blue (uj c > uja) detuning. 

FIG. 2 (a) S-wave scattering cross section an ( in cm 2 ) 
as a function of energy in MHz or momentum fiP\ (in 
1CP 22 gm cm/sec) for two Cs-atoms sharing an optical 
excitation in a cavity without losses (r c = 0, Tr = 2 
MHz) and with loss (r c = 5 MHz, dashed) with param- 
eters as in Fig.l. (b) Idem, for Rydberg Cs-atoms shar- 
ing an excitation near to a = 600GHz, ka — IbOKHz, 
kb = 0.99ka, and u c — loa = 1-OKHz for an ideal cav- 
ity (solid line) and dissipative cavity with r c = 2KHz 
(dashed lines). 

FIG. 3 Spectrum (arbitrary scale) of | Xi^2v) —*\ 
9A9B, 1} transition for X and n symmetry of the quasi- 
binding state with T e ff — 8MHz and randomly oriented 
diatomic axis. The parameters are as in Fig.l. 
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